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Abstract —This  paper  presents  results  of  several  campaigns  at  the  Aredbo  Observatory  using  the  new  HF  high 
power  ‘heater’  fa^ity.  The  specific  measurement  pertains  to  detecting  VLF/ELF  radiation  from  the 
ionospheric  dynamo  current  system  that  has  been  periodically  modified  by  pulsing  the  HF  high  power 
transmitter  at  VLF  or  ELF.  B^use  of  the  nonlinear  relationship  between  ionoqrhnic  conductivities  and 
electron  temperature,  the  high  power  HF  source  can  produce  modification  of  the  ionospheric  conductivities 
and  consequently  of  the  dynamo  current  system.  That  this  eSKt  indeed  does  occur  is  verified  by  detecting  the 
VLF/ELF  radiation  from  this  current  source.  Most  of  the  data  presented  in  this  paper  are  concerned  with 
measurements  made  at  ground  level  beneath  the  modified  current  region.  Limited  data  are  presented  for  ELF 
signals  received  in  Pennsylvaiiia,  U.S.A.,  over  a  long  path  from  the  HF  *heater’  facility  at  Tromse,  Norway. 


I.  INTRODUCTION 

All  elementary  description  of  the  generation  of  ELF 
and  VLF  radiations  by  HF  modification  of  the 
ionospheric  current  system  is  first  presented.  This 
phenomenon  involves  the  irradiation  of  the  iono¬ 
spheric  plasma  by  a  strong  HF  wave  which  is 
modulated  at  some  frequency,  fl  The  periodic  {riasma 
heating  that  results  from  this  wave  is  responsible  for  the 
emission  of  radio  waves  at  a  frequency  (1  from  the 
ionospheric  current  system;  hence  this  effect  is 
commonly  referred  tp  as  nonlinear  demodulation  or 
ionospheric  detection. 

The  basic  physics  of  the  nonlinear  demodulation 
effect  is  iUusf|iNRee^Fip  1  andlThedectrons,  within 
the  area  of  plasma  irradiated  by  the  HF  beam, 
experience  periodic  heating  at  the  HF  modulation 
frequency.  Knee  the  ionospheric  conductivity  is 
elertroa  temperature  depen^t,  the  conductivity 
undergoes  a  similar  pmiodic  variation.  Natural 
ionospheric  currents(e.g.  dynamo,  polar  eiectrojet,  and 

Cuatorial  dectro^)  which  pass  through  the  heated 
lion  are  modulat^  by  the  conductivity  changes. 
iJader  ambioit  conditions  the  deneity  of  the  iono¬ 


spheric  current  is  given  by 

J©  =  eo®o«  (i) 

where  Eq  is  the  strength  of  the  geoelectric  field  which 
drives  the  currents  and  oo  =  is  the  ambient 
conductivity  tensor.  Within  the  heated  region,  the 
current  density  is 

J-efT’JEo'Jo+MT’JEo.  (2) 


Ai»(7’J»e(r.)-eo-  (3) 

Thus,  HF  modulation  of  the  ionospheric  conductivity 
pves  rise  to  a  periodic  current  at  the  modulation 
frequency,  Q,  of  density 

Jo-A^rjEo-  (4) 

Thn  current  fr  a  source  of  G-fiequency  radiation, 
radiated  by  a  novd  type  of  antenna  syston. 

Nonlincw  demodulation  phenomena  of  the  type 
described  above  were  obseived  during  the  Soviet 
Iwating  experiments  at  Oorkii  (Obthantsev  et  oL, 
1974).  During  these  experiments  the  pedar  curmts 
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Fig.  1.  Geometry  of  the  nonlinear  demodulation  experiment. 


were  modulated  at  frequencies  from  1  to  7  kHz  using 
a  5.75  MHz  carrier  with  an  effective  radiated  power  of 
IS  MW.  Signals  were  regularly  observed  during 
the  daytime  at  2.5  and  4  kHz  but  not  at  1  or  7  kHz. 
The  field  strength  at  2.S  kHz  varied  from  2  x  10'^  to 
2S  X  10'^  /t\  m~‘.  A  distinct  diurnal  variation  of  the 
signal  was  reported,  with  the  intensity  reaching  a 
maximum  at  noon  and  being  generally  unobserv¬ 
able  at  nighttime.  Subsequent  measurements  by 
Budiun  et  ai.  (1977)  have  shown  that  the  average 
height  of  the  region  of  generation  of  the  nonlinear 
current,  Jq,  is  about  75  km. 

The  above  results  were  for  magnetically  quiet 
conditions.  Getmantsev  et  al.  (1974)  also  report  that 
measurements  performed  under  magnetically  dis¬ 
turbed  conditions  show  that  the  VLF  field  is  highly 
variable  during  these  conditions  and  that  the  signal 
strength  sometimes  greatly  exceeds  the  intensity  of 
signals  received  during  magnetically  quiet  conditions. 
Nighttime  measurements  by  Kapustin  et  al.  (1977) 
during  geomagnetic  disturbances  show  that  the 
observed  VLF  signal  strength  is  well  correlated  with  the 
strength  and  position  of  auroral  currents  as  detected  by 
groundbased  magnetometers.  The  VLF  signals  were 
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Fig  1  Diagram  showing  ELF/VLF  modulated  ionoapheric 
currant  by  HF  heating 


absent  during  magnetically  quiet  nighttime  conditions. 

Recently,  Turunen  et  al.  (1980)  and  Stubbe  (1980) 
have  reported  the  observation  of  VLF/ELF  signals  due 
to  nonlinear  demodulation  of  MF  and  HF  waves 
during  geomagnetic  disturbances  in  northern 
Scandinavia.  Stubbe  (1980)  observed  a  maximum 
signal  strength  near  2.5  kHz. 

X  THEOKY  OF  VIjr/ELF  RADUTION 
2.1.  Heating  theory 

If  one  neglects  ‘'.i  effects  of  electron  thermal 
conductivity  and  plasma  transport,  the  electron  energy 
balance  equation  can  be  written  simply  as 


where  7^  is  the  electron  temperature,  n,  the  electron 
density,  Kb  the  Boltzmann’s  constant,  P,  the  rate  of 
electron  thermal  energy  production  due  to  HF  heating 
and  (dl//dt),j  is  the  rate  of  electron  thermal  energy  loss 
due  to  the /th  collisional  process  between  electrons  and 
neutral  particles.  It  is  assumed  that  the  electron  and 
neutral  particle  energy  distributions  are  Maxwellian 
and  characterized  by  temperatures  7;  and  T„ 
respectively,  with  T,  =  T,  prior  to  the  introduction  of 
HF  heating 

When  a  region  of  the  ionosphere  is  iHufflinated  by  an 
HF  source  radiating  at  frequency/,  the  local  powm  ^ 
density  5  at  height  h  is  given  by 

S(h)  -  So(^J  exp[^-2*£y(fr0  dA' 

where  So  P-iG/4id^  the  power  dmuity  at  some 
refermoe  bei^t  Ao  below  the  base  of  the  kmotpheie, 
and  Pt  is  the  transmittn’  power  (lectanpilar  hMtiiig^ 
poise  assumed),  0  the  antenna  gain  over  isotrt^'^: 
radistion,  A  the  free  qwce  propagation  ocmstant,  and ' 
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X(h')  the  absorption  index  at  height  h'.  The  HF  source  is 
assumed  to  be  radiating  vertically  in  a  narrow  beam 
over  which  the  antenna  gain  is  essentially  constant. 

Differentiating  equation  (6)  with  respect  to  height, 
one  obtains 

(7) 

The  first  term  on  the  right  accounts  for  the  decrease  in 
power  density  as  the  HF  beam  spreads  with  increasing 
height,  while  the  second  term  represents  the  power 
dissipation  per  unit  volume  as  a  result  of  absorption  by 

,  electrons.  One  can,  therefore,  write 

t  P.  =  2kxS.  (8) 

As  a  heating  pulse  propagates  through  the 
ionosphere,  energy  is  absorbed  from  it  and  converted  to 
electron  heating  as  dictated  by  equation  (8).  Electron 
heating  modifies  the  absorbing  properties  in  the  region 
of  propagation,  since  x  is  a  function  of  7^  through  its 
dependence  on  the  electron-neutral  collision  frequency 
v(7l).  Thus  the  heating  pulse  causes  additional 
absorption  on  itself. 

Electron  heating  results  for  a  strong  field  have  been 
given  in  great  detail  by  Tomko  (1981).  These  results 
were  obtained  numerically  from  equations  (5),  (6)  and 
(8);  a  few  salient  features  can  be  stated  here.  After 
sustained  CW  heating,  electron  temperatures  can  be 
increased  by  about  700  K  at  70  km  during  a  daytime  D- 
region,  with  an  effective  radiated  power  (ERP)  of  100 
MW  usinga  5.0  MHz,  x-mode  HF  signal.  Under  pulsed 
conditions,  at  70  km  due  to  a  250  ps,  5  MHz  x-mode 
signal,  the  electron  temperature  reaches  a  steady  state 
in  less  than  250  ps  and  the  plasma  cools  somewhat 
faster  than  it  heats  up  for  100  MW  ERP.  The  available 
power  density  at  70  km  remains  constant  with  time. 
This  is  not  the  case  in  the  upper  D-region  (90  km) ;  here 
the  power  density  drops  abruptly,  from  its  initial  value, 
to  a  steady-state  value  on  a  time  scale  of  about  500  ps. 
This  effect  is  due  to  self-absorption.  The  electron 
temperature  at  90  km,  on  the  other  hand,  responds 
much  more  slowly  than  the  power  density,  requiring 
about  3.5  ms  to  reach  a  steady-state  level.  Self¬ 
absorption  severely  limits  the  available  heating  power 
in  the  upper  f>-regiott. 

2.2.  Ionospheric  conductivities 

The  periodic  temperature  variations  due  to  pulsing 
the  HF  heater  at  a  VLF/ELF  rate  can  alter  the 
ionospheric  conductivities  in  two  ways.  On  short-time 
scales,  the  collision  frequency  is  modified,  while  on 
longer  timescales(i.e.  >  10~*s)tbedectrondeittity  will 
be  modified.  The  numner  in  which  S' depends  upon 
and  n„  the  dectron-neutral  collision  frequmey  and  the 


Ills 


electron  density,  is  given  by 


MF.)  =  n.(r.) 


v«,(j;) 

0Jg  +  vL(7;)’ 


=  n.(T.) 


e^  cuq 
m.  o>g  +  v^(T.)’ 


=  n.(rj 


1 

w.  v„(T,)’ 


(9a) 

(9b) 

(9c) 


where  Op,  Oh,  and  o^  are  the  Pedersen  (or  Transverse), 
Hall,  and  Longitudinal  conductivities,  respectively.  In 
this  work  it  is  assumed  that  the  heating  periods  are 
short  enough  that  electron  density  modifications  do 
not  occur. 

The  relationships  between  electron-neutral  collision 
frequency  (s'  *)  and  electron  temperature  (K)  as  given 
by  Banks  and  Kockarts  (1973)  arc 

vn,  =  9.32  X  10'  '"iin,  (1  -3.44  x  10-’T.)T„  (10a) 

Vo,  =  1.21  X  10'  ‘°noj  (1+2.15  x  lO'^T.^'^lT.*'*, 


Vo  =  5.49xlO'‘®noT,*« 


(10b) 

(10c) 


where  Hn,,  Ro,  and  Rq  are  the  density  of  N2,  Oj  and  O, 
respectively  (cm  '  *). 

Because  of  the  limited  vertical  extent  of  the 
ionospheric  £-region,  charges  accumulate  at  the 
boundaries  of  the  conducting  layers ;  thus  the  resulting 
vertical  component  of  dynamo  current  will  cease  to 
flow.  The  current  then  is  horizontally  directed.  Under 
these  conditions  the  conductivity  tensor  can  be  reduced 
from  a  3  X  3  tensor  to  a  2  x  2  tensor.  Using  x  and  y 
coordinates  for  magnetic  south  and  east,  the  layer 
conductivity  can  be  expressed,  following  Rishbeth  and 
Garriott  (1969),  as 


where 


and 


ffi  sin*  J+op  cos*  /’ 

ULffii  sin  I 
Ol  sin*/ + Up  cos*/’ 


(12a) 

(12b) 


mth  /  being  the  magnetic  dip  angle. 

Thus  periodic  temperature  variations  will  cause 
paiodk  variations  in  0^,,  and  ip  turn 
modulating  the  horizontal  currait  syston. 
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2.3.  VLF/ELF  raditaion  concept 

Estimates  of  the  radiated  fields  observed  at  ground 
level  directly  underneath  the  modulated  currents 
follow  rather  directly  from  classical  antenna  theory. 
The  spatial  volume  of  modified  current  fills  the  region 
of  plasma  that  is  intersected  by  the  main  beam  of  the  HF 
heater  antenna  pattern.  This  is  approximated  as  shown 
in  Fig.  3  by  an  ‘inverted’  pyramid.  Heights  ff  b  tuid 
define  the  vertical  extent  of  the  modified  currents  which 
in  turn  depend  upon  the  values  of  a„  <r^  and  o,,  under 
ionospheric  and  heating  parameters,  llie  region  from 
H B  to  Hi  is  decomposed  into  current  sheets  each  having 
a  thickness  of  1.0  km.  The  approximate  E  and  H  fields 
calculated  for  ground  level  can  be  easily  expressed 
(Richardson,  1982)  and  summed  over  all  current 
sheets  between  and  Hi.  The  radiation  computer 
code  thus  will  allow  numerical  evaluation  of  these  fields 
once  the  modified  currents  as  a  function  of  altitude  are 
computed. 

2.4.  Computational  methods 

Sections  (2.1)-(2.3)  briefly  described  the  necessary 
computational  steps  required  in  order  to  estimate  the 
signal  intensities  of  VLF/ELF  radiation.  In  our  work, 
the  entire  process  was  modeled  numerically  using  three 
computer  codes : 

(a)  chemical  model  computer  code. 


(b)  electron  heating  computer  code, 

(c)  radiation  computer  code. 

Other  work,  by  Stubbe  and  Kopka  (1977),  for  instance, 
developed  analytical  expressions  for  the  VLF/ELF 
radiation  by  a-wuming  that  temperature  variations  can 
be  expressed  as  e^.  In  the  present  approach  numerical 
methods  are  employed  in  order  to  iiwlude  self¬ 
absorption  effects  and  temperature-dependent  time 
constants. 

The  basic  steps  taken  to  calculate  the  ELF/VLF 
fields  are  as  follows :  from  the  10.7  cm  solar  flux  and 
solar  X-ray  data  on  a  given  day  for  the  latitude  and 
longitude  of  the  heater  site,  the  chemical  model  code 
creates  the  electron  density  distribution  for  the  D-  and 
lower  £-regions  (see  Rowe  et  al.,  1974).  The  output 
from  this  code  is  used  in  the  dectron  heating  code 
(Tomko,  1981),  along  with  parameters  of  the  heating 
facility  sudi  as  fi'equency,  polarization  mode,  radiated 
power,  pulse  width  and  repetition  rate,  to  create  a  file  of 
time  history  of  electron  temperature  as  a  function  of 
altitude.  Also  calculated  are  the  values  of  a„,  and 
(r,y  vs  time  and  altitude.  From  these  values  the 
horizontal  currents  are  calculated  on  the  basis  of  a 
dynamo  electric  field  of  2S  mV  m~  *  directed  in  the 
east-west  direction.  Finally  the  radiation  code 
(Richardson,  1982)  calculates  the  expected  total  field 
at  ground  level.  In  this  calculation  the  spatial  volume  of 
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modulated  current  is  assumed  to  be  in  free  space.  The 
free  space  field  is  multiplied  by  an  absorption  term  for 
the  pertinent  VLF/ELF  signal  to  approximate  the 
actual  field. 

X  EQUIPMENT 

3.1.  Description  of  VLF/ELF  receiving  facility 

A  block  diagram  of  the  receiving  system  used  for  this 
work  is  shown  in  Fig.  4.  The  receiving  system  is  a  tuned 
ELF  receiver  employing  synchronous  detection.  The 
in-phase  and  quadrature  outputs  of  the  synchronous 
detector,  along  with  the  A.G.C.  voltage  from  the  heater 
monitor  receiver,  are  sampled  by  the  data  acquisition 
unit  and  recorded  on  magnetic  disk.  The  A.G.C.  voltage 
of  the  heater  monitor  receiver  is  recorded  in  order  to 
define  the  heater  on-off  status. 

A  balanced  loop  antenna  1  m  in  diameter  consisting 
of 200 turns  of  No.  1 4  wire  and  having  a  tuned  Q  of  64  at 
2.S  kHz  was  used.  A  0.5  m  diameter  one  turn  calibration 
loop  was  located  coaxially  at  a  distance  of  0.5  m  from 
the  main  loop.  The  current  through  the  calibration  loop 
was  measured  with  the  aid  of  a  50  ohm  resistor  (riaoed  in 
series  with  the  loop.  The  calibration  loop  was  used  to 
generate  a  known  field  intensity  for  calibration  of  the 
receiving  system. 

A  Princeton  AppUed  Research  model  113,  low  noise 
high  input  impedance  (100  Mohms)  differential  input 
preamplifier  was  used  in  the  system.  The  gain  of  this 


preamplifier  is  normally  fixed  at  40  dB  with  a  3  dB 
bandwidth  of  1-10  kHz  when  measurements  are  made 
at  2.5  or  5  kHz.  The  second  filter  and  amplifier  in  the 
system  has  a  bandwidth  of  only  +  2%  /„  where  f  is  the 
center  frequency,  and  a  rolloff  rate  of  1 1 5  dB  per  octave. 
The  gain  of  this  amplifier  is  adjustaUe  from  zero  to  60 
dB  in  order  to  compensate  for  variations  in  the  received 
signal  strength.  A  Princeton  Applied  Research  model 
5204  synchronous  detector  with  the  input  sensitivity 
adjusted  to  2.5  mV  full  scale  was  used.  The  reference 
channel  input  to  the  detector  is  derived  from  a 
frequency  standard  accurate  to  two  parts  in  10’;  this 
stability  allows  the  accurate  measurement  of  phase 
changes  in  the  ionospheric  signal  and  also  allows 
integration  times  of  up  to  30  s  to  be  used  at  the  output  of 
the  synchronous  detector.  A  time  constant  of  30  s, 
corresponding  to  a  3  dB  bandwidth  of  0.005  Hz,  is 
normally  used  in  data  recording. 

Prior  to  the  July  1981  observations,  the  receiver  was 
not  accurately  calibrated,  but  for  tbe  July  and 
December  1981  data  the  receiver  was  carefully 
caliinvted.  Calibration  of  the  receiver  is  based  on  the 
following  relationships.  Jean  et  al.  (1961)  give  the 
mutual  impedance  between  two  coaxial  one  turn  loops 


c[^(l-y)K(k)-E(k)j,  (13) 


Fig.  4.  Block  dfagraai  of  the  ELF/VLF  nesiving  sysiHL 
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where 


4r,r,  T'* 

+r.)^+4*J  ’ 


r ,  and  rj  are  the  transmitting  and  receiving  loop  radius, 
respectively,  d  is  the  distance  between  loops,  w  is  the 
angular  frequency  of  operation  and  K{k)  and  E(k)  are 
elliptic  integrals.  The  r.m.s.  voltage  Vj  produced  across 
the  receiver  loop  by  a  current  7,  in  the  calibration  loop 
will  then  be 


V2  =  f.ZuA^e.  (14) 


where  N  is  the  number  of  turns  and  Q  is  the  tuned  Q  of 
the  receiving  loop.  Using  Faraday’s  law,  the  voltage 
which  will  appear  across  the  receiving  loop  due  to  a 
traveling  electromagnetic  wave  will  be 


yiT  =  -jQcoBAN  cos  <l>,  (15) 


where  at  is  the  angular  frequency  of  the  wave,  B  is  the 
r.m.s.  magnetic  flux  density  (Wb  m~^),  ^  is  the  angle 
between  B  and  the  axis  of  the  loop,  and  A  is  the  eflective 
area  of  the  loop  (m^).  Taking  the  magnitude  of  ^'2T■ 
assuming  =  0  and  setting  Kj  =  1^2t<  one  obtains 


|fl|  = 


wA 


(16) 


Thus  one  obtains  the  relationship  between  the 
magnitude  of  the  current  in  the  calibration  loop  and  the 
magnitude  of  a  traveling  electromagnetic  wave  which 
produces  the  same  voltage  across  the  receiving  loop. 
Since  B  =  ftH  equation  (16)  becomes 


|f/|  = 


|f»|Z.2 

coAfi 


(17) 


where  |7f|  is  the  equivalent  magnetic  field  intensity 
(A  m~‘).  Upon  substituting  the  proper  constants 
into  equations  (13)  and  (IS)  one  finds  that  for  the 
antenna  equation  (17)  reduces  to 


|H|  =  !7, 1  X  8.6x10-*.  (18) 


The  actual  calibration  of  the  system  proceeds  as 
follows ;  a  large,  approximately  1  V,  signal  is  placed 
across  the  calibration  loop  and  series  resistor.  The  level 
of  this  signal  is  accurately  measured  and  the  current 
through  the  loop  calculated.  The  calibration  signal  is 
then  attenuated  by  80  dB  and  the  equivalent  magnetic 
field  intensity  calculated  from  equation  (18).  A 
recording  several  minutes  long  is  then  made  of  the 
known  calibration  signal.  This  procedure  is  repeated 
for  each  frequency  of  operation  and  at  calibration 
signal  levels  other  than  80  dB  below  the  1  V  reference 
level  in  order  to  check  for  any  frequency  dependence  in 
the  calibration  or  nonlinearity  in  the  receiver. 

During  the  analysis  of  ELF  data  the  amplitude  of  the 
rMeived  data,  with  appropriate  corrections  for  changes 
in  receiver  gain,  is  compared  with  the  amplitude  of  the 
calibration  recording  in  order  to  ascertain  the  field 
intensity  of  the  ELF  signal. 

3.2.  Parameters  of  the  heating  facilities 

The  facilities  at  Arecibo,  Puerto  Rico  and  Tromso, 
Norway  have  been  used  to  mcxlulate  the  ionospheric 
currmit  system.  Most  of  the  data  presented  in  Section  4 
of  this  paper  were  obtained  at  the  Arecibo  Observatory. 
A  few  long-path  measurements  made  at  University 
Park,  Pennsylvania,  relied  upon  the  Tromse  facility. 
Table  1  summarizes  the  important  parameters  asso¬ 
ciated  with  each  ionospheric  heating  facility. 

4.  OBSERVATIONAL  RESULTS 
4.1.  Measurements  at  the  Arecibo  Observatory 

VLF/ELF  data  were  observed  over  four  heating 
campaigns  at  the  Arecibo  facility.  Meastirements  were 
made  directly  beneath  the  modified  current  region  at 
daytime  and  nighttime.  Table  2  summarizes  important 
information  about  the  campaigns,  such  as  dates,  HF 
heater  parameters,  and  received  signals.  The  last  two 
columns  of  Table  2  give  an  indication  of  the  solar  and 
geomagnetic  activity  at  times  of  VLF/ELF  reception. 


Table  1.  Ccmparison  of  the  Arecibo  and  Tromse  HF  ionospheric  heating  focilities. 


Arecibo,  Puerto  Rico 

Tromse,  Norway 

Geographic  location 

1«’’N,67”W 

TIFN,  1»°E 

Magnetic  latitude 

3rN 

6rN 

Magnetic  dip 

51" 

TT 

Frequency 

3-12  MHz 

Z75-8MHZ 

Power 

800  kW 

1.2-1.4  MW 

4  X  8  array  of  orthogonal 

6x6arrayof 

non-piaaar  log-pwiodic 

orthogonal  Out 

dipole  arrays 

dipoleB 

Antenna  zenith  gain 

23-26  dBi 

24  dBi 

Antenna  beamwidth 

y-up 

8.5" 

Eflhctivc  radhtted  power 

160-320  MW 

300-350  MW 
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Table  2.  Summary  of  data  and  parameters  for  four  campaigns  at  the  Aredbo  Observatory 


Campaign 

1  Dates 

Heater  frequency 
(MHz) 

Heater  power 
(MW  ERP) 

VLF/ELF 

frequencies 

(Hz) 

Solar  10.7  cm  flux 
(10-”Wm-*Hz-‘) 

Range  of 
hourly 
equatorial 

Dst  valims 
(gamma) 

1 

1-3  December  1980 

3.175 

75 

4000,5000 

197.9-213.3 

-16  to  -40 

2 

24-30  March  1981 

3.175, 8.175 

150,84 

5000,2500, 

191.4-208.2 

8  to  -56 

1000,500 

3 

20-23  July  1981 

3.170, 3.175 

100, 125 

Very  weak 

signals. 

228.1-235.9 

6to  -15 

received 

sporadic^y 

4 

30  November- 

3.170 

75, 100 

500^  2500, 

230.3-292.4 

6  to  -23 

13  December  1981 

99Z032, 

574.712 

The  rest  of  this  section  will  discuss  in  more  detail  the 
characteristics  of  the  received  VLF/ELF  signal. 

A  typical  data  set  contains  in-phase  and  quadrature 
voltage  signals  recorded  using  the  receiver  system 
described  earlier.  These  signals  are  then  converted  into 
amplitude  and  phase  before  being  processed  for 
possible  indication  of  VLF/ELF  signals.  At  the  present, 
the  data  are  processed  using  the  fast  Fourier  transform 
algorithms  of  Beroland  and  Dolan  (1979).  While  the 
VLF/ELF  waves  are  generated  by  pulsing  the  HF  hi^ 
power  transmitter  at  a  VLF/ELF  rate,  the  transmitter 
is  turned  on  and  off  periodically  in  order  to  facilitate 
further  the  detection  of  weak  signals.  The  on-off  periods 
are  typically  either  4  or  10  min  in  duration. 

In  applying  FFT  data  processing,  special  care  was 
exercised  in  optimally  choosing  the  sample  length  in 
view  of  the  extremely  weak  signals  which  were  being 
received.  If  the  sample  length  is  limited  to  one  on-off 
period,  a  momentary  noise  burst  could  lead  to  a  false 
indication  of  the  fundamental  component.  On  the  other 
hand,  if  it  contains  too  many  on-off  periods,  it  oouM 
obscure  the  short-term  presence  of  true  ELF  signals. 
Thus  the  data  were  processed  with  a  sample  length  of 
three  on-off  pniods  for  data  taken  with  the  4  min 
period ;  a  sample  length  of  two  on-off  periods  was  used 
for  data  takoi  with  a  10  min  period. 

Figures  shows  a  sample  ofprooessed  data  obtained  . 
by  pulse-modulating  the  3.170  MHz  heater  HF  carrier 
at  a  frequency  of  2.5  kHz  with  a  4  min  on-off  period ;  the 
heater  was  on  for  2  min  and  then  off  for  2  mia.  Figures 
5(a)  and  (b)  show  the  ampbtude  and  the  phase  of  the 
s^g^  respectively,  whereas  Fig.  5(c)  shows  the  spectral 
conqwnents  obtained  by  applying  the  FFT  analysis  to 
the  ampBtnde  of  the  signal.  The  times  when  the  healer 
was  on  are  indicated  alo^  the  thariMse  of  Figs  5(a)  and 
(b)^  and  commence  at  2020^2(04  and  2028111.  The  2.5 


kHz  amplitude  clearly  alternates  between  on  and  off 
periods  indicating  the  presence  of  2.5  kHz  signals.  It 
should  also  be  noted  that  the  phase  signal  is  well  locked 
at  a  nearly  constant  level  during  the  heater-on  period, 
while  it  is  fluctuating  randomly  during  the  heater-off 
period  due  to  absence  of  a  signal.  This  feature 
demonstrates  the  capability  of  the  measurement 


Fig.  S.  Amplitude,  phase  (rad)  and  spectrwn  of  aaqilitude 
fliMuatkmsrff  2.3  kHz  ngaidsoMasuiad  at  Afocibo  from  2017 
to  2028  UT  on  3  Deoeadisr  1981. 
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technique  to  observe  the  phase  variation  of  VLF/ELF 
signals. 

Presence  of  a  2.S  kHz  radiated  signal  is  also  evident 
from  Fig.  5(c)  where  the  fundmnental  frequency 
component  of  4  min  on-off  period,  4.17  mHz,  and  its 
harmonic  components,  are  shown  with  thicker  lines. 
The  fundamental  frequency  component  is  dominant ;  it 
is  interesting  to  note  that  odd  harmonic  components 
prevail,  a  feature  which  is  explained  by  the  near-square- 
wave  time  functions.  Additional  data  are  shown  later  in 
Fig.  7,  describing  Talse’  tests,  in  which  the  presence  of  a 
strong  ELF  signal  at  2.S  kHz  is  also  evident. 
Furthermore,  it  should  be  noted  that  these  data  shown 
in  Fig.  7  were  observed  continuously  over  an  hour  near 
noon  whereas  the  data  shown  in  Fig.  S  were  recorded 
during  late  afremoon. 

4.2.  Long-path  measurements 

In  January  198 1  an  attempt  was  made  to  detect  long- 
path  ELF  signals  at  the  Pennsylvania  State  University, 
University  Park,  Pennsylvania.  The  origin  of  the  ELF 
generation  was  Tromso,  Norway,  where  Dr  Peter 
Stubbe  of  the  Max-Planck-Institut  fur  Aeronomie  was 
investigating  local  ELF  phenomena  due  to  polar 
electrojet  modulation.  This  effort  took  place,  over  five 
days,  with  3-8  h  of  data  each  time.  Since  it  was  a  first 
attempt  at  recording  long-path  signals,  the  measure¬ 
ments  encompassed  several  different  modes  of 
operation:  on-off  mode  at  fixed  ELF  as  described 
earlier  in  Section  4.1,  on-off  mode  with  ELF  stepped 
periodically  over  three  different  frequencies,  CW  mode, 
etc. 

Prelimiiury  results  obtained  for  an  8  h  test  on  11 
January  1981,  are  presented  in  this  section.  The  main 
objective  of  the  measurements  was  to  evaluate  the 
feasibility  of  long-path  detection  of  ELF  waves,  using 
an  on-off  mode  with  a  10  min  period. 

The  FFT  analysis  shown  in  Fig.  6(a)  depicts  the  time 
variation,  throughout  an  8  h  period,  of  tte  strength  of 
the  frmdamental  frequency  component,  1.67  mHz,  of 
the  10  min  on-off  period.  Figures  6(b)  and  (c)  show  the 
spectral  distribution  of  two  data  samples  taken  at  21 14 
s^  1932  UT,  respectively.  The  fundamental  com¬ 
ponents  corresponding  to  these  two  samples  are  also 
identified  in  Fig.  6(a).  Figures  6(b)  and  (c)  stow  not  only 
a  distinct  fundamental  frequency  component  of  1.67 
mHz  but  also  strong  odd  hanmmic  components  as 
were  seen  in  Fig.  5(c).  These  two  results  are  presented 
not  only  because  they  diow  distinct  ELF  signal 
components  but  also  because  the  times  m  whidi  they 
were  taken  are  beheved  to  coincide  with  enhanced 
poiif  slectrojet  activity  seen  at  Tromsn. 

Ahhoogh  at  present  only  a  portion  of  the  kmg-patb 


u,i  a'50 

pi 

K  n  l/W 


0.000)  l|  »|  I  |l  )  l|  I,  I  ,1  ,1  (  ■! 

0  Z  4  6  8  O  12 

FREQUENCY  (mHi) 


FREQUENCY  (mHt) 

Fig.  6.  Temporal  variation  of  (a)  amplitude  of  fundamental 
(XHnponent  and  (b)  and  (c)  spectrum  of  amplitude  fluctuations 
demonstrating  the  on-off  period  for  2.073  kHz  signals  received 
from  Tromse  during  the  period  1505-2336  UT  on  1 1  January 
1981. 


data  have  been  processed  the  results  obtained  so  far  are 
encouraging. 

4.3.  ‘False’  test  investigations 

Numerous  possible  sources  of  frdsdy  generated 
'lignals  have  been  investigated.  Two  of  the  most 
probable  sources  which  were  extensively  studied  are : 
(l)reception  of  ELF  sidebands  radiated  by  the  heating 
transmitter  and  (2)  nonlinear  demodulation  of  the 
beating  ngnal  in  the  ELF  receiver. 

An  examine  of  the  possible  generation  of  &lse  signals 
caused  by  the  sidebands  of  the  heating  transmitto'  is 
seen  by  choosing  the  modulation  frequency  to  be  IS 
kHz  when  the  HF  beater  frequency  is  3.170  MHz;  in 
this  case  the  heating  trannnitter  could  possibly 
produce  a  sideband  at  IS  kHz.  Altboi^  the  order 
(n  -  1267)  of  thn  particular  sideband  is  quite  huge, 
the  heating  transmitter  is  extremdy  powerfrd  and  the 
ELF  reoriver  is  very  snative,  iJlowing  for  the 
posnUe  reception  of  radiated  skMwnds. 

In  order  to  test  for  this  effect,  Um  following 
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experiment  was  performed  several  times.  The  heating 
transmitter  was  operated  at  3.170  MHz  with  2.S  kHz 
square-wave  modulation  until  a  reasonably  strong 
ELF  signal  was  positively  received  asshovra  in  Fig.  7(a). 
Then  the  frequency  of  the  heating  transmitter  was 
changed  to  3. 171  MHz  and  the  signal  shown  in  Fig.  7(b) 
was  observed.  At  this  new  heating  frequency  all  possible 
sidebands  ]noduced  by  the  heater  are  well  outside  the 
bandpass  of  the  synchronous  detector  (0.005  Hz).  Since 
Figs  7(a)  and  (b)  show  essentially  no  change  in  the 
strengfo  of  the  received  ELF  signal  for  3.170  and  3.171 
MHz,  it  is  concluded  that  no  false  effect  due  to  the 
radiation  of  sidebands  from  the  heating  transmitter  is 
contaminating  the  actual  received  signab. 

A  second  source  of  false  signals  in  the  system  could  be 
nonlinear  demodulation  of  the  heating  signal  in  the 
ELF  receiver.  In  order  to  establish  that  the  system  is 
rather  immune  to  this  type  of  effect  three  steps  were 
taken.  First  the  threshold  level  of  HF  signal  strength 
required  at  the  input  of  the  preamplifier  to  produce  a 
detectable  false  effect  was  determined ;  second  the  level 
of  HF  signal  present  at  the  input  of  the  preamplifier 
during  an  actual  experiment  was  measured  and,  third, 
filter  traps  at  the  heater  frequency  were  installed  into 
the  system.  The  HF  traps  used  were  made  removable  so 


Fig.  7.  Amplitudes  of  ZS  kHz  signals  for  studying  iUse  signal 
dbcu  at  Aiedbo  during  the  period  ISOO-IMO  UT  on  9 
Dember  1981. 


that  the  ELF  signal  could  be  observed  both  with  and 
without  the  traps  in  place. 

In  order  to  determine  the  level  of  HF  signal  required 
at  the  preamplifier  input  to  produce  a  false  signal  at  the 
output  of  the  detector,  a  balanced  signal  source  at  the 
heater  frequency  was  injected  into  the  input 
preamplifier  of  the  ELF  receiver  with  a  50%  duty  cycle 
.and  a  pulse  repetition  frequency  of  2.5  kHz;  it  was 
determined  that  6.3  mV  peak-peak  or  larger  must  be 
present  at  the  input  in  order  to  produce  a  detectable 
false  effect  with  the  gain  of  the  system  set  as  for  normal 
sign^  recording. 

During  the  July  1981  observations  the  actual  HF 
signal  level  at  the  input  to  the  preamplifier  was 
measured  and  found  to  be  of  the  order  of  1.0  mV  peak- 
peak  maximum.  Thus  the  HF  signal  received  by  the 
loop  antenna  should  not  cause  any  significant 
nonlinear  demodulation  effects  to  occur  in  the  ELF 
receiver. 

Although  the  tests  outlined  above  would  reasonably 
preclude  any  nonlinear  demodulation  in  the  ELF 
receiver  an  additional  test  was  made  in  which  HF  traps 
tuned  to  the  heater  frequency  were  connected  across  the 
inputs  to  the  preamplifier.  The  traps  which  were  used 
consisted  of  open  circuited  coaxial  transmission  lines 
one-quarter  wavelength  long  at  the  heater  operating 
frequency  and  produced  a  measured  HF  attenuation  of 
30  dB. 

Figures  7(b)  and  (c)  show  the  results  of  ELF  signals 
received  without  and  with  the  traps  connected 
respectively.  The  small  variation  in  the  signal  levels 
depicted  in  Figs  7(b)  and  (c)  is  due  to  a  slow  ionospheric 
fade  and  to  a  small  change  in  the  tuned  characteristics  of 
the  receiving  antenna  after  installation  of  the  filters. 
Hence  it  is  concluded  that  there  was  no  significant 
demodulation  of  the  transmitted  HF  signal  by  the  ELF 
receiver. 

$.  DISCIJSSION 

This  paper  has  shown  that  it  is  possible  to  detect  the 
ELF/VLF  radiation  that  results  from  modulating  the 
dynamo  current  system  by  HF  heating  of  the  Z>-  and 
lower  £-layers  of  the  ionosphere.  Several  campaigns  of 
observations  at  the  Aredbo  Observatory  have 
rqieatedly  shown  that  radiation  from  500  to  5000  Hz 
can  be  recdved.  This  paper  has  examined  the  pos«bility 
of  receiving  falsely  created  signals,  and  condudes  that 
the  signals  received  are  indeed  generated  within  the 
ionosphere.  A  few  long-path  (Norway-Pennsylvania) 
propagation  results  have  been  described;  here 
modulation  of  the  polar  electrt^  was  induced  from 
Troina0,Norway,daring  experiments  bdngcoDdiicted 
by  the  Max-Planck-Institut  ffrr  Aeronomie.  The  FFT 
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analyses  indicate  with  reasonaUe  certainty  that  the 
ELF  radiation  propagated  over  this  long  path. 
Considerable  amounts  of  data  were  recorded  during 
those  experiments  and  are  undergoing  further  analysis 
to  evaluate  the  statistics  associated  with  this  successful 
detection. 

The  field  strengths  which  were  received  cannot, 
unfortunately,  be  conqrietely  verified  yet  by  the 
antenna  model  described  in  Section  2;  this  model 
predicts  field  strengths  several  times  those  that  were 
actually  measured.  It  is  believed  that  predicting  the 
modified  currents  from  a  change  in  conductivity  is  not 


accurate,  and  that  the  inductance  the  current  path 
may  need  to  be  included. 
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